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Monterey, California. The purpose was to visually analvze the flowfield around the 
model while simulating ship motion in rough seas, with a correspondingly modelled 
Open-ocean atmospheric boundary laver. 

The two flow visualization techniques utilized were neutral density helium 
bubbles and a liquid aerosol, each technique providing varving degrees of on-body and 
off-body analysis. Limited success was achieved with the aerosol, but the helium 
bubbles produced excellent photographic results and an accurate visualization of the 
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I. INTRODUCTION 


Helicopter flight operations aboard UL. S. Navy combatants and auxiliary ships 
have become an important aspect of present fleet readiness. These operations can 
become quite hazardous, particularly in the takeoff and landing phases when combined 
with excessive ship motion and superstructure airwake turbulence. The Naval Air Test 
Center (NATC) presently determines the safe operating envelopes for various ship- 
helicopter combinations, and according to Carico, McCallum, and Higman [Ref 1] 
there is a backlog involving eleven different helicopters and twenty classes of ships. At 
this time, disregarding the inevitable development of new ships and helicopters, it 1s 
estimated that all the safe operating envelope combinations cannot be determined 
mami the remainder of this century. Thus, it is highly desireable to develop a 
simulation technique accurately establishing these safe envelopes, which would save 
time, costs (now at $75,00U to $150,000 per combination), and require no operational 
assets. 

This report is concerned with the analysis of the ship airwake turbulence of an 
oscillating model by flow visualization with respect to the above problem. Two flow 
visualization techniques Were attempted to accomplish this: neutral density helium 
bubbles and an aerosol spray, each of which provided various levels of insight towards 
the flowfield. Each method also required different photographic and lghting 
procedures to achieve optimum results. Experimentation was conducted in the Naval 
Postgraduate School's low speed wind tunnel modified to simulate an open-ocean 
atmospheric boundary laver velocity profile and associated low altitude turbulence 
levels. An oscillating mechanism was designed and modified by Prof. J. V. Healey, and 
was used to approximate ship motion in relatively rough seas. The ship model used 
Was a generic destroyer constructed of a medium-density polystyrene foam body with 
an alumunum bow, scaled approximately at 170:1 to a Spruance-class destrover. Both 
on-body and off-body flows were investigated by the flow visualization techniques 
listed above. Flows were investigated initially in the single degree-of-freedom modes of 
roll and pitch, then were analyzed in combinations of the above to possibly determine 


anv non-linear coupled affects. 


This report deals with a portion of an ongoing investigation. Preceding research 
included verification of the atmospheric boundary layer modelled in the wind tunnel, as 
Well as evaluation of a sinular generic model in static conditions only, using various 
flow visualization techniques. It is intended that further study will analyze the flows 
about specific ship types, using highly detailed models (some of these include the 
Spruance-class destroyer, Enterprise-class carrier, and New Jersev-class battleship). 
The final phase of research anticipates the inclusion of specific helicopter simulators to 
study the interactions of the two bodies in the dynamic open-air ocean environment. 

Photographic results are included and discussed in the final chapter. as Well as 


reconimendations for the follow-on phases of the project. 
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Il. ATMOSPHERIC BOUNDARY LAYER SIMULATION 


A. THEORETICAL CONSIDERATIONS 

Numerous studies have indicated that the turbulent boundary laver simulation 
Paced solciveon the mean velocity profile is insufficient. According to Healey 
[Ref. 2: p. 46], there are four important parameters which must be considered: 


1. The mean longitudinal windspeed. U, which 1s the instantaneous windspeed 
averaged over a specific length of time. 


2. The standard deviation. o, of the windspeed fluctuations about the mean which 
vields the “turbulence intensity” when divided by the mean windspeed. 


Lia 


The longitudinal scale length of the turbulence, L, a measure of the mean length 
Gmuleeamest ellergetic eddies in the turbulence. 


Seeeihe Spectrum {Unction of the turbulence, an indication of the frequency 
distributions of the turbulence energy. 


These four parameters have been empirically shown to be functions of the mean 
windspeed, the elevation in the boundary laver, and the roughness length scale, the 
latter being a surface irregularity coefficient not directly related to ground obstacle 
height. The roughness scale is designated as Z¢). measured in meters. and ranges from 
10> meters for a glassy surface to 3-4 meters for verv rough-textured terrains. A more 
@erailed listing for various types of terrain is available in [Ref. 3] 


Garratt (Ref. 4] approximates the roughness length Zp as 
10 exp {-0.41 (0.75 + 0.067 x U(10))! ° x 10E-3} (eqn 2.1) 


Where .4]1 1s the Von Karman constant, the square root value represents the neutral 
drag coefficient, and U(10) designates the mean windspeed at a 10 meter altitude. 
Keep in mind that this is a general approximation, since several parameters mav be lost 
due to experimental scatter. The 10 meter elevation figure is used in the equation 
because this is a standard reference height often used by meteorologists, and it also 
coincides with the mean height for the helicopter landing decks of modern frigates and 
destrovers. 

Recent experimentation has revealed that the turbulence intensity has some 


dependence on the mean windspeed, but the effect 1s small enough that it will be 


neglected here. \fore importantly, it 1s the fluctuations in the longitudinal velocity 
component that determine the turbulence. The two major physical parameters 
influencing these velocity fluctuations are the altitude and roughness length scale. As 
discussed above, the altitude will be considered constant at 10 meters and a mean 
windspeed of 25 msec (about 50 kts) is indicative of “normal” operating conditions in 
heavier seas. 

Using the mean windspeed value of 25 m sec in equation 2.1, Zg is found to be 
0024 meters. Comparing to other sources [Refs. 3,5], the turbulence intensity 1s 
calculated to be 0.14. The term “rough seas” here is used in the same context as that 
defined by Kent [Ref. 6: p. 22], with a mean wave height of I4 feet. For this condition 
(sea state 5), the roughness length will be considered to range from .001 to .01 and the 
corresponding turbulence intensity range to be 0.13 to 0.17. Recall from earlier, the 
turbulence intensity 1s the standard deviation of the fluctuations divided by the mean 
wind speed. 

More important than the longitudinal turbulence length scale parameter is the 
ratio of this length scale to the ship's beam. As this ratio increases, the actual flowfield 
about the ship appears to be less random and more time-dependent. 

Several spectrum function models are currently available to describe the 
atmospheric turbulence, such as those bv Kaimal, Dryden, and Von Karman. The Von 


Karman function 1s more popular and 1s given by 
i ye Ree 56 . 
S(n) = (4Lo) ; U{1+ 70.8(n }7} (eqn 222i 


Where S(n) 1s the value of the Von Karman spectrum function. n’ 1s defined as Ln,U, L 
is the turbulence length scale, U 1s the longitudinal mean velocity, and n is the 


turbulence frequency. 


B. EXPERIMENTAL SIMULATION 

Numerous approaches have been developed in atmospheric boundary-laver 
simulation (or ABL), but few take into account all four of the parameters listed 
previously. The best source of information available in this area has been found to be 
the Journal of Wind Engineering and Industrial Aerodynanucs. One major difficulty in 
simulating an accurate over-ocean ABL 1s in the comparison with experimentally 
recorded data. since the ABL 1s different for changes 1n wind velocity, sea state, and 


meteorological conditions. 


Also, the source of data being compared must be examined. NATC data 1s 
generally obtained using ship anemometers as the primary source of windspeed 
measurement. These devices measure only the relative windspeed (between the actual 
windspeed and the ship’s speed), and as stated earlier, the atmospheric boundary laver 
changes only with respect to the true windspeed. Of less importance, but a justifiable 
consideration of error source, is that NATC tests have shown anemomieter positioning 
in ship superstructures 1s suscepuble to airwake interference and therefore the reliability 
of the instruments themselves is questionable. Thev are normally calibrated to avoid 
this, but meteorologists would not use them 1n a data collection study. 

In dealing with ship motion relative to the ABL, scaling factors must be utlized 
to truly model the combined airwake in the wind tunnel. It is here the Strouhal 
number is introduced. This is a dimensionless quantity defined by the frequency of 
motion oscillation, n, multiplied by the characteristic length, L, (in this case, the ship's 
beam is used) divided by the speed, L, as given in the following equation: 


~~ 


Sipouhal Number = nxLl/tl (eqn 2.3) 


While n, L, and LU may vary from the actual system to the wind tunnel model, their 
ratios, defined bv Strouhal numbers, must be equal for accurate modelling. Typically, 
as was the case in this situation, the model characteristic length is fixed bv the size of 
the test section, and L may often be liunited bv the nunimum Revnolds number needed 
and by the power capabilities of the motor driving the fan. This leaves the model 
oscillation frequency as the primary variable and thus is easily determmned. The 
calculations used for the Strouhal Numbers are tabulated in the following chapter in 
Table I. 

In view of the above, the most successful technique encountered in developing a 
reliable ABL wind tunnel model was presented by Counihan [Ref. 7]. He utilized 
vortex generators well forward of the test section to create the appropriate mean 
Windspeed profile. Following the vortex generators Were roughness elements, mounted 
low to the wind tunnel floor, to provide the necessary low-level turbulence. The mean 


speed profile, as described by Arva [Ref. 8], is the well known logarithmic relationship 


ee — (1) K)in( Zz) (eqn 2.4) 
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where U’ ts the frictional velocity, K is Von Karman’s constant (previouslv used in 2.1 
as Q.41), Z is the elevation and Z, is the roughness parameter (.001 to .O1, as previously 


determined). This expression may be refined into an exponential form as 
re a — / ies 
Ua (zen Z,) (eqne2ss} 


where now the roughness parameter becomes the gradient height, Z, in this form. For 
the selected mean speed of 25 m sec at 10 meters elevation for this experiment, the 


mean Velocity profile in 2.5 becomes 

L = 25x(Z; 10)! (eqn 2.6) 
Through the correct combination and placement of the vortex generators and the 
roughness elements, the above velocity profile with verified turbulence intensities 


(representing the atmospheric shear stresses) was obtained in the Naval Postgraduate 


School's low-speed wind tunnel [Ref. 9: pp. 21-28]. 
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Hil. SHIP MOTION 


A. THEORETICAL CONSIDERATIONS 
Actual ship motion is a complex, non-linear system with six degrees of freedom 
(or DOF). The system may be approximated by the sixth-order equation of the 


following mass-damper svsteim: 


M Y(t) + DY(t) + K Y(t) = F(t) (eqn 3.1) 


where M, D, and K are six-by-six matrices and F and Y are six-dimensional vectors of 
(x.v,z,9,0,y). It is noted that F(t) may be random or sinusoidal forcing functions in 
Mier lor Of fOrceS OF MOmMeENtS. (X.v,z) represent the three translational modes of 
surge, sway, and heave (respectively), while (@,0,w) describe the roll, pitch, and vaw 
rotational modes. The system is made more complex in that coupling between modes 
occurs; normally roll, swav and vaw are coupled together as are heave, pitch and surge. 
Finally, ship interaction with the immediately surrounding water and wave action 
affects the mass (1), damping (D). and restoring (K) matrices so that the final svstem 
requires verv large-scale computation and numerical approximation methods. Perhaps 
Heeemost difficult calculation is referred to as the eneounter frequency.” This 1s a 
function of the ships natural motion frequencies (one in each DOF) and the 
frequencies of the ocean waves, each having their own spectrum functions: it also has 
a strong dependence on the angle between the ship’s heading and the dominant wave 
direction. 

Two distinct methods currently prevail in ship motion prediction programs. The 
more classical tume-domain approach calculates (or approximates) the coefficient 
matrices, and often utilizes coordinate transformations to decouple certain modes for 
more simple analyses. These are often required to use small-angle approximation 
theory and thus the accuracy falls off as higher amplitudes of motion are encountered. 
The frequency analysis method correlates the energies of the encounter wave to the 
ship's response. the resultant scaling factor referred to as the Response Amplitude 
Operator, or RAO. Determination of these operators is relatively simple, but a huge 
number is required since the RAO changes with ship's heading and speed, sea 
condition, wave frequency, and each DOF. 


ie) 


Healey [Ref. 2: p. 52] lists three ship-motion programs in use, the one from David 
Tavlor Naval Ship Research and Developement Center (DINSRDC) being the most 
sophisticated. The following assumptions are applicable: 

e Monohull ship form (entails all present U.S. Naval Ships) 


e Strip theory (a common aerodynamic technique analyzing the ship's transverse 
sections in two dimensional flow) 


e Linearized equations are valid (an update to the program accounts for non- 
linear roll effects) 


e The RAO’s are derivable from wave and ship motion harmonics 
e Heave, pitch, and surge are uncoupled from roll, sway and vaw. 

This is a relatively effective ship motion program, but adjustments must be made 
as more information in this field becomes available, particularly in the area of non- 
linear angular motion. It is also imperative that this, as with any simulation 
techniques, be validated as soon as time permits against the actual ships which are 


being modelled. 


B. EXPERIMENTAL SIMULATION 

It is the intent of this phase of the project to establish an oscillating mechanism 
for a generic ship model, thus several approximations and assumptions were made 
initially. Primary emphasis was on the individual analysis of the major modes of heave, 
pitch, and roll and any interactions between those. This is supported by Kent 
[Ref. 6: p. 41], in which he states each of the six motions are analyzed only one at a 
time, then the attempt 1s made to trace the effects of each mode on all the others. Due 
to fluid resistance forces, even for a ship in calm waters, any movement in any one 
mode will cause a response in the ship in all five other modes. Simplifying assumptions 
were made for each mode, such as the following for the pitching mode [Ref. 6: pp. 
- 76}: 


(1) The swell consists of a regular series of waves, uniform in dimension and 
sinusoidal in shape 


AN 


7 


if 


(2) The “encounter angle” is zero degrees (or head-on) at constant speed 
(3) Pitching is isochronous 

(4) Pitch damping from water resistance is neglected 

(3) Ship wake and ocean wave interference is neglected 


(6) Variations in virtual mass and effects from ship’s rise or fall is ignored. 
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The other primary objective at this stage is to optimize the flow visualization and 
photographic techniques for this oscillating, bluff body-tvpe flow. With these goals in 
mind, an extensive ship’s motion program as described previously is not required at 
this time, since a more basic understanding of the aerodynamics involved is being 
sought. This provides for a simple ships’ motion model indeed, but much more 
research is required before accurate ship motion programs can be extensively applied to 
wind tunnel testing. 

As discussed earlier on page 13, rotational speeds were determined by equating 
the Strouhal Numbers, shown by equation 2.3, for the ship and the generic model. The 
simulated conditions previously listed are summarized here in Table 1, as well as the 
parameters used in the wind tunnel. (Ship characteristic length and roll rates were 


obtained from personal conversation with Mr. J. Higman of NATC). 


Wea pio lees 
PREQOUENGIES OF MOTION 


MODEL 


17 meters Q.1 meter 


L beam 


Velocity 2 iisec 2.0 mysec 
Roll Rate 10-1 58sec; cycle 

equates to 0.067-0.1 Hz vields 0.9-1.4 Hz 
Pitch Rate 7-10 secrcvcle 


equates to 0.1-0.14 Hz wields 1-4-1,9 11z 





Bluff bodies, as referred to above, are defined as non-streamlined bodies 
generating separated flow over a substantial proportion of their surface. These bodies 


are often characterized by sharp leading-edged geometric shapes of wide angles. Until 
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recently little concern has been shown in their relation to the surrounding flowfields. 
This area of study is now referred to as “bluff bodv aerodynamics,” and the primarv 
interests lies in the fields of architecture and automotive engineering. The model used 
during this testing may be considered sinular to a long, low building, thus those 
governing principles may be applied to the project. Unfortunately, for anything other 
than a basic geometric shape like a cylinder or a rectangle, little has been accomplished 
in this area with respect to detailed analvtical testing. Bearman has done extensive 
study on an oscillating cvlinder, but concedes that more work must be done. 
[Ret 1G7 pple. a 

Two phenomena particularly interesting with oscillating bluff bodies is the 
generation of trailing vortices and flow detachment. It would be advantageous in 
certain situations, such as a helicopter landing on a ship, to have information like the 
amplitudes of the trailing vortices, areas more prone to detaching flows and therefore 
more turbulent flow, and how these areas and their frequencies of occurrence are 


affected by bluff body oscillations, for some exaniples. 
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IV. EXPERIMENTAL APPARATUS 


A. WIND TUNNEL FACILITIES 

The low-speed smoke tunnel located at the Naval Postgraduate School in 
Monterey, California, was used for all experimentauon. It consists of a 15 square foot 
intake section with a honevcomb screen combination to prevent debris ingestion and 
provide some directionality to the inbound flow and reducing the turbulence in the test 
secuon. A square bell-tvpe contraction cone leads to a 5 square foot, 22-foot long 
closed test section, partially enclosed by plexiglass sheets in the walls and ceiling for 
observation and photography. A constant-speed, AC-motor 1s used to drive a variable 
pitch 5-foot diameter fan which is preceeded by protection louvres, located directly 
behind the test section. The sixteen-blade fan 1s capable of producing speeds from 
virtually O to approximately 35-feet per second. The tunnel exit section is then routed 
to the environment for exhaust. 

Several modifications were required to obtain a proper ABL simulation. The 
scheme suggested by Coumbhan (et al) was implemented in the following manner: 
mortex eenerators, to slow the lower-level flow and create the desired velocity profile. 
were constructed and installed at the leading edge of the tunnel. Four were found to 
be sufficient in producing the profile (as per verification measurements described in the 
following section), but spaces between each were wide enough to allow “tubes” of 
unaffected flow, which required the addition of the three tapered delay cones. The 
desired turbulence intensity was then obtained by using three-eighth inch dowels 
ranging from 2 to 8 inches in height as the roughness elements. Equally spaced one 
inch blocks were attempted earlier with poor results. The dowels were randonuy placed 
and the turbulence intensity was measured, then several iterations were required before 


the desired intensity level, as described in Chapter II, was achieved. (Fig. 4.1). 


B. VERIFICATION OF RESULTS 

Bolinger (Ref. 9: pp. 20-25] completely describes the verificauon procedures and 
methods, which are summarized here. Measurements were obtained from two different 
ietruinents tO enstire accurate data was taken: first by an EDM-2500C 
micromanometer and then by a Dantec hotwire anemometer attached to a single wire 


probe. One hundred forty-four data positions were measured in the test section area, 
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probe. One hundred forty-four data positions were measured in the test section area. 
with sixteen horizontal and nine vertical coordinates. It is noted here that this thesis 
deals primarily with flow visualization techniques and its associated theory. A 
concurrent report written by LCDR Gary Selman deals with the data acquisition 
portion, thus any further details in this aspect may be obtained from that source. 


iaeres 2 and 3 {rom Reference 9 are included here for continuity. 


C. OSCILLATING MECHANISM 

[It was decided to simulate only the highest-amphtude motions, as was discussed 
earlier, separately and then to investigate anv non-linear coupling effects between 
modes. Thus the mechanism was designed to sustain heave, pitch, and roll or 
combinations of heave with roll and pitch with roll. The device incorporates two 
“Scotch Yoke” mechanisms connected by two joined shafts, and two motors (one for 
rolling motion. the second driving the yokes for pitch and or heave motions), which are 
mounted to a frame and centered just below the test section. 

A five-amp, | 2-Hp. Minarik motor controlled the vertical motion, using a belt to 
Rotate tO COniiected shafts, in turn driving the two yokes mounted to either end of the 
Medel. lhe Scotch yoke consists of a cylindrical flange with an eccentric pin, 
adjustable along the radius of the flange. The pin slides along the channel of the voke, 
mhich is driven vertically by the pin as the flange 1s rotated (see Fig. 4.2). The vokes 
Were maintained in the true vertical bv Teflon guides. The two shafts connecting the 
vokes were connected at a coupling to facilitate transition from heave to the pitch 
modes. They were easily disconnected at the coupling and could be rotated through 
90-degree phase increments, where in-phase represented pure heave and 1$0-degrees 
delineated the pure pitch. A pitch and heave combination mode was available at the 
90- and 270-degree phase positions but was not analvzed for this experiment. 

The roll mode was driven by a 1 $-horsepower, |-amp Bodine motor, again using 
an eccentric cylindrical flange, but here it was connected directly to a control arm 
which fed into a pin through the center of the model. Both the flange and the arm 
length could be adjusted to provide desired roll angles, ranging from approximately 5- 
to 20- degrees. The roll arm connection to the model was by a simple pin, which was 
easilv removed to shift to non-rolling modes. 

As alluded to earlier, the model was to be analyzed in a flow with variety of vaw 
angles, ranging from zero degrees (head-on to the bow) to 35 degrees either side of the 


bow (a full range of motion through 360 degrees was available, but does not practically 
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simulate typical at-sea operations). It was desired that the entire oscillating mechanism 
could be rotated, and so it was mounted on a columin attached to the center of the 
frame, with a locking handle to prevent any undesired rotation during the experimental 


Pus. 


D. THE SHIP MODEL 

In this preliminary phase of analvsis, it was decided to use a generic destrover 
model with relatively simple geometric shapes to simulate the deck levels and the 
superstructure, vice a more specific model type such as those listed in the introduction. 
In this way it was hoped to gain some foresight in coupled-mode effects and general 
characteristics of trailing vortex generation, particularly directly above the after-end 
area Where most helicopter platforms are installed. 

The following desireable characteristics Were considered in deternuning the model 
composition: 


(1) ease of construction (since it was originally envisioned to use numerous 
generic models, examining affects of subtle changes of geometry) 


(2) durability (to withstand a variety of oscillation frequencies up to three Hertz) 


(3) hghtweight (for minimal strain on the mechanism’s motors and to reduce 
effects of dynamic loading) 


(4) low cost (always a driving factor) 

Also, it was desired to vield easily to modifications, such as the addition of vent holes 
for smoke aerosol injection. Of the materials considered, polyurethane foam of 
medium density (8 Ibs cubic ft) met all the characteristics listed. Surfaces sanded 
relatively well with a sealant applied, and a flat black lacquer could then be directly 
spraved to the exterior. This low-reflective coating was desired for optimum contrast 
during photography, and was therefore also used inside the entire test section area of 
the wind tunnel. 

The platform consisted simply of a 5-foot diameter plywood disk cut out of the 
tunnel floor at the center of the test section. The model shape was then cut into the 
disk with approximately a half inch gap for model clearance in the rolling mode. A 
minor obstacle in implementing this svstem occurred when attempting to prevent 
anibient air from entering the test section through the gap between the model and the 
cutout in the disk. Commercial tvpes of bristle and closed-cell polvurethane 
weatherstripping proved inefficient due to lack of compressibility. The conunercial 


foam with a self-adhesive strip also tended to peel off the cutout edges during extended 
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periods of oscillation. The best results were obtained with a lower density (open-cell) 
polyurethane (the type found in seat cushions) attached to 1-1 4-inch wide batting 
strips along the cutout edge. A good contact cement was used, and the batting 
provided a larger surface area for bonding than did the 3.4-inch unmodified disk. This 
provided a sufficient barrier to unwanted flow caused by the pressure differential during 
the wind tunnel operations. Finally, the disk was reinstalled and supported by circular 
sections of steel plate mounted to the underside of the tunnel floor, making it easy to 
rotate the model disk oscillating mechanism svstem through any desired vaw angles. 
The disk was also clamped at three positions from underneath to ensure it was 
maintained level with the remainder of the test section floor. 

Both the motors were controlled by variable-speed Muinarik voltage control 
boxes, which were niounted on the outside of the tunnel wall but inside the observation 
deck (see Fig. 4.1). The desired rotational speeds were verified by a Monarch 
Instruments TACH IV electronic tachometer, which operates be directing a collimated 
light beam onto a reflective marker attached to either the yoke (for measuring the pitch 
mode) or the outside of the cvlindrical flange (for roll). The reflected pulses are 
compared against a 6-\{Hz. internal crystal tumebase, and the resultant speed 
ineasurement is displaved directly in RPNf on an LED readout. The instrument 
incorporates a microprocessor to compute the time interval between successive 
reflected pulses, which were received be either directly pointing the instrument at the 
reflective markers or by using a remote sensor. The last option was the most preferred 
since it Was not required to leave the observation room to take the measurements 
under the wind tunnel, and the controllers could then be adjusted with real-time 
feedback. 
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V. EXPERIMENTAL METHODS 


A. FLOW VISUALIZATION TECHNIQUES 
Numerous techniques are available to the aerodynamicist for the enhancement of 
qualitative flow study. As with all engineering situations, each has its own benefits and 
drawbacks depending on the requirements of the problem. Flow visualization may be 
approached from two major areas, that of on-body analvsis and that of off-body, or 
the area surrounding the bodv in study. A common dilemma to both areas of on-bodv 
and off-body analvsis is that it 1s imperative that a negligible disturbance to the 
ii@eeeld be introduced to ensure accttrate results are obtained. As an example, 
compare the flow visualization aspect to the data acquisition portion of the project. A 
pitot-static probe may be used nearly anywhere in the test section (except very close to 
the model and tunnel walls due to mutual interference) to determine the freestream 
merocit\, Dileany area downstream of the probe cannot be evaluated concurrently 
Mess the exact disturbance effects fromthe probe are known and can be corrected for 
in the data reduction. Likewise, any flow visualization techniques introduced to the 
flow. whether on-body or off-bodyv, must not alter that which it 1s attempting to 
analvze. 
1. On-body Flow Techniques 
This is the area immediately surrounding the model, and all current techniques 
utilize some sort of application. Techniques considered here included oil and clay 
smears and tufting. The smears would provide a long term record of the local 
streamlines, photography would be less critical and the flow around the entire model 
could be analyzed in detail, but at the low speeds used in this study the required 
streaking would probably not occur. Previous work by LCDR Bolinger for the static 
model entailed the utilization of very thin fluorescent nylon filaments (0.7-thousanths 
of an inch in diameter). These tufts were attached in an evenlv spaced in a grid along 
the model, with application procedures detailed in [Ref. 9: pp. 32-33]. Due to their 
small size and flexibilitv they had virtually no influence on the flowfield, and the 
dvnamics of the local streamlines were easily traced, as seen in Fig. 5.1. Due to delavs 
In inodel construction and implementation, this technique was unable to be used with 


the oscillating model. 





Figure 5.1 Example of Fluorescent Minitufts, Static Model. 


‘“e 


2. Off-body Flow Techniques 
_ These methods typically employ the injection of light-scattering particles into 
the flow, acting as tracers to the flow streamlines. The particles should be of a 
sufficiently small size, at most around a hundred microns, so as to introduce no 
disturbances to the flow but must still be able to disperse the light for visualization 
purposes. Iwo methods were used. 

An exception to the particle size was employed, consisting of neutrally- 
bouyant helium bubbles ranging in sizes of one to three millimeters in diameter. 
Mueller [Ref. 11] provides a detailed description of a similar bubble generating system 
and lists several applications he encountered during his research in flow visualization 
methods. This technique was attractive because bubbles could trace individual 
streamlines of the flow, and also because the quick dispersion of the substance in 
moderate turbulence was not experienced as it is in smoke techniques. Helium bubbles 
were created using a Sage Action, Inc. (SAI) system, which consisted of a bubble 
generator console, an ejector head, and a vortex generating chamber (see Fig. 5.2). 
The console controlled the bubbles size, density and generation rate by a mixture of 


soap-film solution and helium, with 50-psi air used to drive the bubbles through the 


28 


KA NA | compressed 
ify NJ Ly N] 


bubble 


helium 
tank 


air source 





Gonerol 
console 





helium, air, 
and BFS lines 
ie ei mixing nozzle 
flexible pate 
tubing 


filter 


chamber 





to nozzle 
inside tunnel 


Figure 5.2 Schematic of Bubble Generating Svstem. 


29 


helium, with 50-psi air used to drive the bubbles through the chamber. The bubbles 
were created at the ejector nozzle and fed directly into the vortex generating chamber, 
where the lighter bubbles impacted the outlet tube in the center of the chamber and the 
heavier bubbles impacted along the walls of the chamber. This allowed the truly 
neutrally-bouyant bubbles to pass through the vortex generator and be directed 
through flexible tubing into the tunnel’s test section approximately 18 inches upstream 
of the model. A nozzle was used to direct bubble flow downstream, and since the pipe 
supplving the nozzle replaced existing dowels in the tunnel it had no effect on the 
turbulence intensity already calculated. The nozzle was supplied from the underside of 
the tunnel, and could be placed at a variety of spanwise locations depending on the 
model vaw angle and the desired area of flow analysis on the model. This system 
produced excellent results for the static model analysis and good results for the 
oscillating model. 

Due to the relatively high Reynold’s Number (about 50,000), another 
technique was attempted in order to obtain more detailed analysis closer to the model, 
particularly the superstructure areas. A smiall particulate substance was desired for the 
finer detail, with the realization that this type of material would be quickly dispersed. 
Smoke was to be avoided due to safety constraints and it was simply messy to work 
with. A liquid aerosol system described by Griffin and Ramberg [Ref. 12] was 
investigated and found to be quite appropriate for the purpose of this project. This 
svstem produced a fine mist which was discharged from strategic areas of the model, 
especially the areas of shed vortices determined from the static analysis. The system 
described by Griflin was easily modified to existing units of the helium bubble 
generator system from SAI, thereby making it easy to implement. The compressed air 
tube was disconnected from the generator console and routed to the inlet of the liquid 
DOP (di(2-ethvlhexyl)-phthalate} agitation chamber; the output was then connected, 
again by flexible tubing, to a jet impactor nozzle inside the drum used as the vortex 
generator for the bubble system. The jet impactor was used to reduce the particle size 
and enable the smaller particles to follow the unaffected flow. The output of the 
chamber then went to a brass manifold with six outlets, from which tubing was run to 
the various areas of the model described earlier. All connections between major 
components of both systems were with flexible tubing except for the aerosol exhaust 
ports which consisted of brass piping through the model and connected to tubing on 


the model’s underside. The DOP was an ideal substance, since it is chemically inert, 


non-toxic, non-corrosive, the system) requires no chemical reaction or combustion. and 


it contains no moving components for very little maintainance. 


B. PHOTOGRAPHIC TECHNIQUES AND LIGHTING 
Each flow visualization technique required somewhat different photographic 
methods and therefore different lighting schemes Were necessary for optimum results. 
With only a few exceptions, the photographic equipment used was a single-lens reflex 
Hasselblad 2000-FCW camera with one of two Hasselblad lenses: either a 110-nim. or 
a 150-mm. The Hasselblad is a professional-format camera, which was desired for 
better resolution on enlargement printing. It became normal procedure to shoot with 
3000-speed, tvpe-107 black and white Polaroid film initially to ensure proper exposure, 
then the final photography was done using Kodak TMAX 400 professional film. To 
account for the difference in film speeds, the TMAX was exposed bv opening the 
aperture one ‘f-stop’ and overdeveloping (“pushing”) the film two stops. Other types 
and speeds of both color and black-and-white films were evaluated, but the best 
contrast Was most consistently obtained with the Kodak TMAX. In all cases. proper 
lighting was the primary obstacle, due tn large part to the restrictive confines of the 
observation photographic deck and the limited size of the plexiglass viewing ports in 
the walls and ceiling. Also, anv wind tunnel is faced with the problems of blockage 
and interference, so discretion had to be maintained in the placement of anv 
illununating devices inside the tunnel itself, particularly upstream of the test section. 
The reference used for all photographic techniques and procedures was Handbook for 
Scientific Photography, by Alfred Blaker [Ref. 13]. The specific techniques emploved 
are now discussed in the following text. 
1. Oscillating Model with Aerosol 
[t was much more difficult to obtain usable photographs when the aerosol 
flow visualization technique was utilized. The lighting procedures first attempted are 
simular to those discussed bv Griflin (Ref. 12: p. 68], in which the camera shutter was 
held open in a blacked-out room and the exposure was determined by the back flash 
set on the opposite side of the tunnel. The flash unit used was a Norman 2000, 
capable of various power levels from: 400 to 2000 W-secs., and the head was flashed 
through a honeycomb screen for directionalitv. The full range of flash power settings 
were also explored. In an attempt to bring out more contrast in the photographs, K-2 
(vellow) and A (red) filters were used for many of the shots, in accordance with Blaker 


(Ref. 13: pp.69-79]. Another effort entailed using Kodak Ectachrome color slide film 


flash techniques combined with fast shutter speeds were initially desired to “freeze” the 
motion of the ship during oscillation, with a wide range of aperture settings and shutter 
speeds examined. 

2. Static/Oscillating Model with Helium Bubbles 

This method produced the best photographic results of the three flow 
visualization techniques used. Since the desired end result was to obtain bubble streaks 
which followed the flow, extended shutter speeds were required. The lower end of the 
scale (fastest shutter speed) was at 1'4-seconds, with any faster speeds than that giving 
short traces and little practicality of analysis. The upper end was at 2-seconds, since 
slower speeds produced so many traces that the overlap was too great to be able to 
follow any one particle. The best results were obtained with a shutter speed of 
l-second for the static case and 1, 2-second for the oscillating runs, with aperture 
settings between F-4 and F-8'11 (depending on the lighting scheme). 

The bubbles showed up best when being illuminated against a totally dark 
background, so the lighting was set up to cast as little light as possible onto the model. 
This became quite challenging at times in view of the fact that the model was three 
dimensional and of irregular shape. The source of lighting which produced the highest 
trace intensities was an EIMN[LAC model R-150-5 high-voltage arc lamp placed inside the 
tunnel and well downstream of the model. This was used to either direct a circular 
beam along the port (left) side of the ship, which in all cases was the 
viewing photographic side, or around the entire model. The latter case worked better 
for the static nose-down runs and some of the oscillation shots. Due to limitations 
with lateral working space inside the tunnel, combined with the arc lamp being housed 
in a three-foot long shroud, the runs at a 35-degree vaw angle required the beam to be 
bounced off of a murror attached to the tunnel wall. Adjustment was time consuming, 
but once in place the effect was substantially improved with the arc lamp. 

The arc lamp, however, provided illumination in limited areas, and interesting 
sections of the flow, such as between the stacks or behind the superstructure, required 
alternate lighting with good directional versatility. It was found that slide projectors 
served this purpose quite well. Two types were used, a Kodal Ectagraphic Model AF-2 
and a 500-W Sawyers “Rotodisc” projector, with three to six being used at anv one 
time depending on the areas of interest surrounding the model. Normally, two or three 
were used on the opposite side of the tunnel to illuminate the areas of the bow and 


stacks, with the remainder positioned at various angles inside the tunnel aft of the test 
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section. An 800-Watt professional theatre lamp was occasionally used, but was bulky 
and was difficult to accurately direct around the model inside the confines of the wind 
tunnel. [he projectors could provide a large area of lighting, sav over the top of the 
model or along the starboard side, otherwise thev were used with slides to create highly 
detailed beams to completely silhouette the ship. These slides were manufactured from 
strips of aluminum tape with the shape of the lighted area being cut out, and since 
construction was relatively easy numerous were made and could be changed at 
convenience (for example, after a change of vaw angle). This method provided 
excellent results for the static cases, but the set-up time was normally quite long and 
laborious. 

The oscillating runs presented more of a dilemma: the lighting could be 
optimized for only one position, but the model motion caused some areas to be unlit at 
Smemearieo! the cycle and producing excessive model glare at other parts. A 
compronuse was made to reduce some of the illumination areas very close to the model 
and accept some occurrences of glare. It is noted here that substantial blur due to the 
ship movement was anticipated but would be accepted since the actual subject matter 
was the bubble traces themselves. Mforeover, it was hoped to determine the ship 


position from the blur through various stages of its cycle. 
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VI. RESORTS 


The flow visualization methods which were utulized in this study have been 
described in detail, along with their associated photographic and lighting techniques. 
The following results are first discussed according to the flow visualization method (as 
in Chapter V), with the helium bubble section subdivided by vaw angles. It is noted, 
and will become obvious as the figures are scanned, that it was the bubble technique 
from which the conclusions were derived. 

The tunnel freestream velocity was held constant at 9.1 feet per second for all of 
the test runs. From the data of Table 2, the freestream velocity can be computed for 
various heights of interest. For example, the velocity over the bow (at z= 2 inches) 
was approximately 6.25 ft’sec., or at the top of the stacks (at z=4.5 inches) it was 
approximately 6.9 ft, sec. 

The pictures presented here are representative of a much larger grouping from 
Which the observations and conclusions were made. These are published solely as an 
example of the results obtained and to help illustrate the validity of the conclusions 
being drawn. It is also noted that much of the resolution from the original 


photographs was lost during the printing process. 


A. AEROSOL INJECTION 

The results for this techmique can only be described as disappointing. In 
attempting to use a high-powered flash from across the tunnel, the results were 
inconsistent and were very sensitive to the angle of the flash head with the honeycomb 
attachment. Without the use of the honeycomb the pictures were extremely 
overexposed, even at the lowest power settings for the flash and the highest aperture 
settings on the camera. It is suspected that the flash head was quite prone to the 
vibration of the tunnel, causing it to become quickly misaligned. 

The system plumbing is also suspected to be a contributor to the frustration 
caused be this method! Eighth-inch tubing was used to transport the aerosol through 
the ports, this may have been too restrictive to obtain a high-density flow. Thicker 
aerosol clouds were easily obtained by increasing the air pressure into the DOP 
agitauon chamber. but this action produced relatively high velocity aerosol jets from 


the exit ports and thus made this trick invalid. Also, the regulator valve for the 
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compressed air inlet was very sensitive to adjustment and a consistent but yet realistic 


aerosol flow was difficult to obtain. 





Figure 6.1 Sample Results with Aerosol. 


The photo included here (Figure 6.1) was taken by the third technique described 
in the photographic procedures section, using Ektachrome color film for higher 
contrast. The difficulty generated from this technique may have been primarily from 
the glare received from the flash bouncing off of the plexiglass viewing port. 

The aerosol injection technique continues to be modified as of the writing of this 
report, including possibilities of installation of larger diameter tubing and shooting the 


photographs from inside the tunnel to avoid any glare problems. 


B. HELIUM BUBBLE STREAKS 
This method of flow visualization became the preferred analvsis technique, since 

once the lighting problems discussed in Chapter V were overcome the procedure was 
straightforward and the results were consistent. At each vaw angle the model was 
photographed at the following positions: 

(1) Zero degrees for roll and pitch (true static condition) 

(ee bullet relly zero pitch. 

(3) Full right roll, zero pitch. 
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(3) Full mght roll, zero pitch. 

(4) Zero roll, full nose down pitch. 

(5) Zero roll, full nose up pitch. 

(6) Oscillating model in both roll and pitch. 
The above positions are defined as follows: full roll (left or right) is approximately 8 
degrees, and full nose (up or down) 1s approximately 4 degrees. All figure titles are in 
reference to the above, in that “nose up” is assumed full nose up at 4 degrees and “left 
roll” refers to full left roll at 8 degrees. The photos of the oscillating model show that 
the ship is blurred, but the helium bubble streaks were generally quite clear. Numerous 
photos show distinct shaded columns above the two stacks. These are from the 
projector lighting at the opposite side of the test section, and were highly pronounced 
only on days when smoke visualization techniques were being used concurrently in 
other nearby facilities. Lastly, little difference could be seen consistently enough in the 
two different roll positions to draw any definite conclusions, thus only one roll position 


for each series of yaw angles is included in the following results sections. 





Figure 6.2 Zero Degrees for Roll and Pitch at Zero Degrees Yaw. 
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Figure 6.4 Nose Up Pitch at Zero Degrees Yaw. 
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Figure 6.5 Nose Down Pitch at Zero Degrees Yaw. 


1. Zero Degrees Yaw 

.The flow generally appears to be quite symmetrical along either side of the 
model, despite the existence of the turbulent boundary laver. Evidence of the 
recirculation zones, although negligible in Fig. 6.2, is seen behind the aft stack in Fig. 
6.3. For the nose up case, not evident in Fig. 6.4, the flow has a very distinct 
detachment over the bow but quickly reattaches and has little affect on the 
superstructure and forward stack. As would be expected, the opposite is true for the 
nose down case (Fig. 6.5), with stronger trailing vortices being generated off the leading 
edge of the superstructure and forward stack. The oscillating shot (Fig. 6.6) has 
substantially fewer bubble traces due to the shorter exposure times, but displays some 
interesting features. A “wave affect” is seen on several traces along the port side, 
presumably due to the motion of the stacks, and the forward stack oscillation is shown 


by the series of inward and outward direction changes of several bubble traces. 
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2. Twenty Degrees Right Yaw 

A foot vortex has generated just aft of the bow and travels along the lower 
port side leading-edge of the superstructure, with numerous weaker vortices trailing 
from various areas of the ship. The initial trailing vortex angles are the same as the 
model's yaw, but as they continue travelling downstream they realign themselves with 
the freestream and join together to form a single large vortex. Figure 6.7 shows the 
shift of the stagnant regions around the front stack. The left roll position (Figs. 6.8) 
shows a stronger trailing vortex at the aft end of the forward stack, again in 
comparison to the zero-degree yaw position, as well as the generation of a trailing 
vortex off of the edge of the superstructure. The pitch positions (Figs. 6.9 and 6.10, 
respectively) display effects similar to those seen at the zero-degree yaw position. The 
oscillating shot (Fig. 6.11) was chosen to show the effect on the recirculation zones: 
bubbles were still captured here and actually maintained oscillations with the model for 


a short time before escaping back to the freestream. 





Figure 6.8 Left Roll at 20-Degrees Right Yaw. 
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Figure 6.10 Nose Down Pitch at 20-Degrees Right Yaw. 
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Figure 6.11 Roll and Pitch Oscillations at 20-Degrees Right Yaw. 
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3. Twenty Degrees Left Yaw 
This position shows primarily a mirror image of the twenty-right yaw position 
(Figs. 6.12 to 6.16). More obvious from this aspect, however, is the development of 
the large corkscrew vortex on the leeward side. Of particular interest are the trapped 
bubbles in Fig. 6.12 directly behind the first stack, and the large area of detaching flow 
off the bow in Fig. 6.14. The bubble traces appear more irregular in Fig. 6.16, 
particularly along the port side, which may be attributable to the ship's oscillations. 
4. Thirty-five Degrees Left Yaw 
Trailing vortices are evident on most of the sharp edges of the model at this 
higher yaw angle (Figs. 6.17 to 6.21). The corkscrew vortex has grown even more in 
size, and the flow generally appears more turbulent throughout the entire flow field. 
The position at right roll (Fig. 6.18) shows the detaching of trailing vortices over the 
entire section before the superstructure. This effect is weakened in the nose down 
position, but much more prevalent with the bow up than has been seen at the previous 
vaw angles. The model in the oscillating mode (Fig. 6.21) displays more irregularities 
in the bubble traces than do the static modes, presumably due to the interaction with 


the ship's motion. 
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Figure 6.15 Nose Down Pitch at 20-Degrees Left Yaw. 
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Figure 6.19 Nose Down Pitch at 35-Degrees Left Yaw. 
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Figure 6.20 Nose Up Pitch at 35-Degrees Left Yaw. 
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Figure 6.21 Oscillating Model at 35-Degrees Left Yaw. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


This report dealt with the investigation of visualizing the flow field about an 


oscillating ship model. As a result of this and an earlier study (Ref. 9), the following 


basic conclusions can be diac 


(1) 


The ship motion simulation was a viable first-step approximation to actual 
ship motion. 


The heltum bubble flow visualization technique was quite successful for both 
the static and oscillating cases. 


The flow freld about both the static model and, to a lesser degree, the 
oscillating model, was effectly analyzed. 


Little determination of coupling effects from the three oscillating modes 
previously listed was realized. This area requires further flow visualization 
study combined with accurate data acquisition techniques. 


The following recommendations are given for consideraton by those who will 


continue with the next phases of this project (as discussed in the Introduction): 


(1) 


ee, 
fad 
Ce 


(4) 


The use of high-resolution video is strongly recommended for a more exact 
study of the oscillaung mode. In addition to enabling a detailed analysis of 
the flow field over a wider range of oscillation frequencies, it would provide 
virtually real-time feedback and bvpass the delays involved in developement 
and printing of photographs. 


The fluorescent minituft technique would provide valuable information in the 
oscillating mode, such as the determination of whether vortex shedding from 


the model body exists. Future models should incorporate tulting as described 
Ineieiees 


Modification to the bubble generator system to possiblv increase the bubble 
output rate would provide the option to photograph at faster shutter speeds. 
This would make the frequency analvsis based on the bubble streaks less 
complicated. 


The author is convinced that the aerosol technique shows enough promuse for 
further modification. One possible avenue, as pointed out in the results 
section, would be to investigate the effects of increasing the tubing diameter 
and use higher air pressures to provide a more dense aerosol. 


The EIMAC high-voltage arc lamp was responsible for the majority of the 
high contrast bubble traces. More illumination sources of this tvpe, with an 
added feature of a geometrically adjustable iris to conform to the model shape, 
ae nected: 
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(3) 


(9) 


(10) 


(11) 


Since lighting was one of the biggest obstacles, alternate illumination schemes 
shoulemeerconsidered. Ihese miaht mclude the use of laser sheets, or LED's 
Or optical fibers from within the ship (this would require no adjustment in the 
static or oscillating case). 


The camera shutter should be svnchronized to the ship oscillations, which 
could easily be incorporated with existing hardware. Thus, if recommendation 
(1.) ts unfeasible at this time (due to funding constraints). this would be an 
alternate method of obtaining a more accurate frequency analvsis. 


‘A switch for a second remote to the electronic tachometer would be helpful to 
facilitate measurement of both modes of motion. 


Some analvsis of the rolling modes may be obscured due to the motion being 
in the same plane as the camera angle. Shooting from the ceiling viewing port 
or even from inside the tunnel mav provide some added insight. The latter 
case would also eliminate problems with reflective glare from the plexiglass. 


Flow visualization is only a first step in accurate flow analvsis, action 1s 
required now to introduce data acquisition methods to achieve more detailed 
Fesmlts it tae models areas of interest. 


Further study is required into incorporating a more realistic ship motion 
program into the NPS low-speed wind tunnel. 
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